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Abstract 

Collider experiments have recently measured the production cross section 
of hard jets separated by a rapidity gap as a function of transverse momen- 
tum and gap size. We show that these measurements reveal the relative 
frequencies for the production of rapidity gaps in quark-quark, quark-gluon 
and gluon-gluon interactions. The results are at variance with the idea that 
the exchange of 2 gluons is a first order approximation for the mechanism 
producing colorless states, i.e. the hard QCD Pomeron. They do qualita- 
tively support the "soft color" or "color evaporation" scheme developed in 
the context of bound-state heavy quark production. 



I. INTRODUCTION: 



Although we do understand strong interactions in the context of QCD, we can at best 
speculate on how to calculate elastic scattering. This is just one example of an interaction 
mediated by the exchange of the "Pomeron" , a state which carries no net color. The routine 
speculation has been that it is, to a first approximation, a state of two colored gluons 
combined into a color singlet. Understanding the Pomeron has been challenging because its 
dynamics is revealed in processes which are not subject to perturbative computation, e.g. 
elastic scattering. The hope has been that it may be instructive to study Pomeron dynamics 
in hard processes that are, at least partially, understood in terms of perturbative QCD: the 
hard Pomeron. Examples include processes involving colorless pairs of heavy quarks, e.g. 
ip^s, or colorless states of light quarks produced in association with a pair of high transverse 
momentum jets: rapidity gaps. In this paper we will show how a treatment of colorless 
states in QCD, suggested by the phenomenology of heavy quark bound states, supports a 
"soft color" model of the Pomeron which is at variance with the idea that it is a structure 
built on a frame composed of two gluons. We will show that recent measurements [|l| of 
the relative frequency for the production of rapidity gaps in quark-quark, quark-gluon and 
gluon-gluon interactions provides qualitative, yet convincing, confirmation of the soft color 
concept. 

The reason why some data on the production of ip- and T-states radically disagree with 
QCD predictions, occasionally by well over one order of magnitude, is that the traditional 
method for performing the perturbative calculation of the cross section is simply wrong 
0. The key mistake is to require that the heavy quark pair forms a color singlet at short 
distances, given that there is an essentially infinite time for soft gluons to readjust the color 
of the cc pair before it appears as an asymptotic ip or, alternatively, DD state. We suspect 
that the same mistake is made in the description of rapidity gaps, i.e. the production of a 
color-neutral quark-antiquark pair, in terms of the exchange of a color neutral gluon pair. 
The ip is, after all, a color neutral cc pair and we have shown p[ in quantitative detail 
that it is indeed produced by exactly the same dynamics as DD pairs; its color happens 
to be bleached by soft final-state interactions. This approach to color is also suggestive of 
the unorthodox prescription for the production of rapidity gaps in deep inelastic scattering, 
proposed by Buchmiiller and Hebecker 0. 

In this paper we emphasize that recent measurements IQ on the formation of rapidity 
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gaps between a pair of high transverse momentum jets shed new hght on the problem of how 
to treat color in semi-hard interactions. When applied to the formation of gaps between 
a pair of high transverse momentum jets in hadron collisions, the "soft color" or "color 
evaporation" approach suggests a formation rate of gaps in gluon-gluon subprocesses which 
is similar, or smaller, than in quark-quark induced events. Consequently, formation of gaps 
should increase with increased transverse momentum, or reduced collision energy of the jets. 
This prediction happens to be antithetical to the one obtained in 2-gluon exchange Pomeron 
models. We show that the data resolve the issue in favor of the "soft color" computational 
scheme and questions the relevance of approximating the exchange of a color-singlet, hard 
Pomeron as a pair of gluons. We also exhibit the predictions of the "soft color" model for 
the production of rapidity gaps at the LHC energy. 

II. ONIUM CALCULATIONS WITH SOFT COLOR: A BRIEF REMINDER 

The conventional treatment of color in perturbative QCD calculations, i.e., the color 
singlet model, has run into serious problems in describing the data on the production of 
charmonium and upsilon states Specific proposals to solve the onium problem agree 
on the basic solution: onium production is a two-step process where a heavy quark pair is 
produced first. At this initial stage all perturbative diagrams are included, whether the cc 
pair forms a color singlet state or not. This is a departure of the textbook approach where 
only diagrams with the charm pair in a color singlet state are selected. In the Bodwin- 
Braaten-Lepage (BBL) formalism [§] the subsequent evolution of the pair into a colorless 
bound state is described by an expansion in powers of the relative velocity of the heavy 
quarks in the onium system. A different approach, the color evaporation or soft color method, 
represents an even more radical departure from the way color singlet states are conventionally 
treated in perturbation theory. Color is, in fact, "ignored". Rather than explicitly imposing 
that the system is in a color singlet state in the short-distance perturbative diagrams, the 
appearance of color singlet asymptotic states depends solely on the outcome of large-distance 
fluctuations of quarks and gluons. In other words, color is controlled by nonperturbative 
interactions. 

In Fig. Ill we show typical diagrams for the production of ^/'-particles representing the 
competing treatments of the color quantum number. In the diagram of Fig. |l|a, the color 
singlet approach, the ip is produced in gluon-gluon interactions in association with a final 
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state gluon, which is required by color conservation. This diagram is related by crossing 
to the hadronic decay ^ ?> gluons. In the color evaporation approach, the color singlet 
property of the ip is ignored at the perturbative stage of the calculation. The ip can, for in- 
stance, be produced to leading order by gg-annihilation into cc, which is the color-equivalent 
of the Drell-Yan process, as shown in Fig. |I|b. This diagram is calculated perturbatively, 
with dynamics dictated by short-distance interactions of range Ax ~ wi^^. It does indeed 
not seem logical to enforce the color singlet property of the ip at short distances, given that 
there is an essentially infinite time for soft gluons to readjust the color of the cc pair before 
it appears as an asymptotic ip or, alternatively, DD state. Alternatively, it is indeed hard to 
imagine that a color singlet state formed at a range m^^, automatically survives to form a 
ip. This formalism represents the original and, as we have shown correct method 
by which perturbative QCD calculations should be performed. 

The evidence is compelling that Nature operates according to the color evaporation 
scheme. The formalism predicts that, up to color and normalization factors, the energy, xp- 
and pT-dependences of the cross section, are identical for the production of onium states and 



DD pairs. This is indeed the case [|TO|,|TT|. Another striking feature is that the production 
of charmonium is dominated by the conversion of a colored gluon into a as in Fig. ^b. 
In the conventional treatment, where the color singlet property of the ip is enforced at 
the perturbative level, 3 gluons (or 2 gluons and a photon) are required to produce a ip. 
Contrary to the usual folklore, ip^s are not produced only by gluon-gluon interactions. As a 
consequence color evaporation predicts an enhanced ip cross section for antiproton beams, 
while the color singlet model predicts roughly equal cross sections for proton and antiproton 
beams. The prediction of an enhanced p yield is obviously correct: antiproton production 
of ip^s exceeds that by protons by a factor 5 close to threshold. This fact has been known 
for some time 10-^. We should note that for sufficiently high energies, gluon initial states 
will eventually dominate because they represent the bulk of soft partons. 

Quantitative tests of color evaporation are made possible by the fact that all ■^/'-production 
data, i.e. photo-, hadroproduction, Z-decay, etc., are described in terms of a single param- 
eter: the parameter determining the frequency by which a charm pair turns into a ip via 
the final state color fluctuations. Once this parameter has been empirically determined for 
one initial state, the cross section is predicted without free parameters for any other. We 



have demonstrated |11] the quantitative precision of the color evaporation scheme by show- 
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ing how it accommodates all measurements, including the high energy Tevatron and HERA 
data, which have represented a considerable challenge for other computational schemes. Its 
parameter- free prediction of the rate for Z-boson decay into ip^s is an order of magnitude 
larger than the color singlet model and consistent with data ||12|| . 

In summary, the soft color approach gives a complete picture of charmonium production 
in hadron-hadron, 7-hadron, and Z decays. The phenomenological success of the soft color 
scheme is impressive and extends to applications to other charmonium and upsilon states 



III. RAPIDITY GAPS AS COLORLESS STATES OF (LIGHT) QUARKS AND 

GLUONS 

We now turn to the implications of the soft color scheme for the dynamics underlying 
the production of rapidity gaps, which refer to regions in phase space where no hadrons 
appear as a result of the production of a color neutral partonic system. The connection 
to charmonium physics is obvious: the ip is a color-neutral cc pair. The important lesson 
from heavy quark phenomenology is that perturbative color octet states fully contribute to 
the asymptotic production of color singlet states, such as ip's. We suspect that this is also 
true for the production of a rapidity gap which represents nothing but the creation of a 
color singlet quark-antiquark pair, as shown in Fig. ^ja for electroproduction. The diagram 
represents the production of final state hadrons which are ordered in rapidity: From top 
to bottom we find the fragments of the intermediate partonic quark-antiquark state and 
those of the target. Buchmiiller and Hebecker proposed that the origin of a rapidity gap 
corresponds to the absence of a color string between photon and proton remnants, i.e. the 
3 X 3 (= 1 + 8) intermediate quark-antiquark state is in a color singlet state. Because color 
is the source of hadrons, only the color octet states yield hadronic asymptotic states. The 
reasonable guess that 

2 -^^g 2 ^ J 

follows from this argument. 

Although this is only a guess, it embodies the essential physics: events with and with- 
out gaps are described by the same short-distance dynamics. Essentially non-perturbative 
final-state interactions dictate the appearance of gaps whose frequency may, possibly, be 
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determined by simple counting. The treatment of color is the same as in the case of heavy 
quark production and leads to similar predictions: the same perturbative mechanisms, i.e. 
gluon exchange, dictates the dynamics of color-singlet gap (tp) and regular deep inelastic 
(open charm) events. 

Our proposal for the (soft) nature of color challenges the orthodox mechanism for pro- 
ducing rapidity gaps sketched in Fig. ^o, where the t-channel exchange of a pair of gluons 
in a color singlet state is the origin of the gap. The color string which connects photon 
and proton remnants in diagrams such as the one in Fig. is absent and no hadrons are 
produced in the rapidity region separating them. The same mechanism predicts rapidity 
gaps between a pair of jets produced in hadronic collisions; see Fig. These have been 
observed and occur with a frequency of order of one percent |1^-|16 . 



The arguments developed in this paper question the hard Pomeron approach: it is as 
meaningless to enforce the color singlet nature of the gluon pair as it is to require that the cc 
pair producing aip is colorless at the perturbative level. Following our color scheme the gaps 
originate from a mere final state color bleaching phenomenon d la Buchmiiller and Hebecker. 
This can be visualized using the diagram shown in Fig. |]b. At short distances it represents a 
conventional perturbative diagram for the production of a pair of jets. Therefore, the same 
short distance dynamics governs events with and without rapidity gaps, as was the case for 
electroproduction. This is consistent with all experimental information. 

Producing a quantitative model for the gap rate may be premature at this point. There 
are at least two consistent interpretations of the present data. The first is based on the 
string picture for the formation of the final state hadrons shown in Fig. ^d. Color in the 
final state is bleached by strings connecting the 3 jet at the top with the 3 spectator di-quark 
at the bottom and vice-versa, resulting in color singlet states at the top and bottom. The 
probability to form a gap can be counted d la Buchmiiller and Hebecker to be 1/(1 + 8)^ 
because it requires the formation of singlets in 2 strings. The data |]T^|-0 is consistent 
with this simple picture which basically predicts that the gap fraction between pp jets is the 
square of that between virtual photon and proton in deep inelastic scattering. One could, 
alternatively, argue that, once color has been bleached between one 3-3 pair, overall color 
conservation will guarantee the color singlet value of the other pair. This leads to a gap 
formation rate which is similar in lepto- and hadroproduction, and can be reconciled with 
the data by introducing a survival probability. The survival probability accounts for the 
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fact that gaps can be filled by the underlying event, e.g. mini-jet production |T^, or by 



higher order processes. The survival probability of the gap is expected to be smaller for 
hadroproduction than electroproduction, thus accommodating the data. 

Discussions of rapidity gap physics have routinely ignored that, besides quark-quark, 
gluon-gluon and gluon-quark subprocesses contribute to jet production in hadron collisions. 
In the gluon-gluon color flow diagram corresponding to Fig. ^ top and bottom protons each 
split into a color octet gluon and a 3-quark remnant in a color octet state. There are now 
(8 X 8)^ color final states. Despite the fact that we can at best guess the non-perturbative 
dynamics, it is clear that the soft color formalism predicts a gap rate which is smaller in 
gluon-gluon interactions. This is in contrast with the 2-gluon exchange diagram of Fig. ^a 
which predicts a gap-rate enhanced by a factor (^|^ in gluon-gluon subprocesses |T^. The 



process is clearly enhanced when replacing the interacting quarks by gluons, because of 
the larger gluon-gluon color coupling. The contrasting predictions can be easily tested by 
enhancing the relative importance of quark-quark subprocesses in the experimental sample. 
This can be achieved by increasing the pt of the jets at fixed energy, or by decreasing the 
collision energy of the hadrons at fixed px- In either case, we anticipate an increased rate for 
the production of gaps in the soft color scheme; a prediction opposite to that of the 2-gluon 
exchange model. We next confront the contrasting predictions with recent data 

Introducing the quantities Fqq, Fqc, and Fgg which represent the frequencies for pro- 
ducing rapidity gaps between a pair of high-p^- jets in hadronic quark-quark, quark-gluon, 
and gluon-gluon collisions, we can write the observed gap fraction as 



where 



da = daQQ + dagc + dace (3) 



represents the decomposition of the cross section for producing laige-Ex jets into quark- 
quark, quark-gluon and gluon-gluon subprocesses. Predictions can be summarized in terms 



of the gap fractions F^j. For the 2-gluon hard Pomeron model |TB| 

9 /9^2 



Fqq : Fqg : Fqc " ^ • 4 • (^4 j • (4) 

In the soft color calculational scheme the Fij are independent of the center-of-mass energy 
and Ft, satisfying 
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Fqq : Fqg : Fgg = a : b : c , (5) 

with c < b < a, in contrast with Eq. (^. Reasonable guesses fall in the range 

(1/9)2 < a < 1/9 , (6) 
ac < b < \fac , (7) 
(1/64)2 < c < 1/64 . (8) 

In order to determine the gap fractions Fjj, we computed the transverse cross sections 
daij/dET in lowest order perturbative QCD, and subsequently fitted the preliminary D0 
data |[l[] shown in Fig. ^. We integrated over the Ft bins given in Table |, imposing that 
both jets have \ri\ > 1.9. The 90% CL bounds on the Fij are 



Fqq - 0.023_o^oii , 



0.010 

11 

0.012 



Fqg = 0.00017^^00 , (9) 



Fgg = 0.0075+°-°°«^ 



-0.0073 



We also show in Fig. ^ the result of our fit. 

The data shows that rapidity gaps are mostly formed in quark-quark collisions with the 
value of Fqq exceeding those of Fgg and Fgq- The 90% CL upper limits on Fqg and Fgg are 
0.012 and 0.016, respectively. The dominance of Fqq simple reflects the experimental fact 
that the fraction of events with gaps increases with Ft as the leading order QCD result for 
{daQQ / dEx) / {da / dET)] see Table |. Clearly, the data is consistent with the soft color model 
for the formation of rapidity gaps and does not support the 2-gluon exchange approximation 
which predicts that processes involving gluons should exhibit larger rapidity gap frequencies. 
The data can be interpreted in terms of the simple color counting previously introduced, 
which predicts Fqq ^ (0.2-0.4)/9 ^ 2-4 10"^ and Fgg - (0.2-0.4)/64 ~ 3-6 lO^^ for a 
survival probability [1^ of 0.2-0.4. This is certainly compatible with our results given the 



uncertain systematics of our procedure. 

Our formalism has non-trivial djTiamics built in. It assumes that the gap fractions Fij 
are independent of any kinematic variables. Therefore, the observed fraction of rapidity gaps 
depends on kinematical variables only through the relative contributions of the 3 subprocess 
cross sections (Tqq, o"qg, and aGG- The observed gap fractions are, for instance, a function 
of Et only at fixed center-of-mass energy, as shown in Fig. |^. The Ft distribution was 
computed using uniform binning at Tevatron and LHC energies. Furthermore, we also 
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predict that the gap frequencies are only a function of the scaled variable xt = Et / y/s. 
Therefore, the gap frequencies are described by a universal curve at all energies; see Fig. 

At LHC energy, the fraction of events exhibiting rapidity gaps associated with moder- 
ate Et jets (< 80 GeV) is quite small (~ 0.6%) due to the large gluon-gluon luminosity. 
Notwithstanding, we anticipate 40 spectacular events for Et > 200 GeV because the gap 
frequency increases to ~ 0.8% in this Et range. On a more practical note, overlapping 
events may make such observations challenging. 

The D0 Collaboration has also measured the dependence of the fraction of events with 
rapidity gaps as a function of the rapidity separation of the leading jets for Et > 30 GeV 



and \ri\ > 1.7 [SQ]. Our predictions are successfully confronted with their recent data |1| in 
Fig. 0. 

The soft color mechanism can also give rise to rapidity gaps with a different ordering in 
rapidity: near the beam are the remnants of an initial hadron, followed by a rapidity gap and 
two hard jets which are separated from the other hadron debris by yet another rapidity gap. 
These events represent diffractive gaps. In the context of the color evaporation model, quark- 
quark and quark-gluon collisions cannot initiate such events because of the impossibility to 
neutralize the color of the 3 (3) remnants with the exchange of soft gluons. Therefore, 
such events only originate in gluon-gluon collisions, and their fraction should decrease with 
increasing Et, or decreasing center-of-mass energy. 

We close with some cautionary comments. Do Figs, ^j-^ suggest that we have formulated 
alternative s- and t-channel pictures to view the same physics? Although they seem at 
first radically different, this may not be the case. Computation of the exchange of a pair 
of colorless gluons in the t-channel is not straightforward and embodies all the unsolved 
mysteries of constructing the "Pomeron" in QCD. In a class of models where the Pomeron 



is constructed out of gluons with a dynamically generated mass [p7f| , pi| , the diagram of 
Fig. ^a is, not surprisingly, dominated by the configuration where one gluon is hard and the 
other soft. The diagram is identical to the standard perturbative diagram except for the 
presence of a soft, long-wavelength gluon whose only role is to bleach color. Its dynamical 
role is minimal, events with gaps are not really different from events without, suggesting 
dynamics similar to color evaporation. Some have argued that in this class of models the 
hard Pomeron is no more than an order a? correction, a view which can be defended on more 



solid theoretical ground ||2^. Others have however challenged the theoretical soundness of 
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this line of thinking [lT9| , p3| . Also note that our discussion is at best indirectly relevant to 
completely non-perturbative phenomena like elastic scattering. There is no short distance 
limit defined by a large scale. The Pomeron exists. 
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TABLE I. Et bins used in Fig. ^ and its respective mean value Et- For comparison, we also 
provide the cross section integrated over the bins for the quark-quark, quark-gluon, and gluon-gluon 
subprocesses. 



Et (GeV) 


Et bin (GeV) 








21.0 


15-25 


60.8 


197 


163 


29.8 


25-30 


3.25 


7.19 


4.02 


35.1 


30-35 


1.12 


2.07 


0.97 


40.3 


35-40 


0.43 


0.68 


0.27 


45.5 


40-45 


0.18 


0.26 


0.084 


50.7 


45-50 


0.079 


0.093 


0.028 


57.4 


50-60 


0.056 


0.057 


0.014 


70.5 


> 60 


0.027 


0.020 


0.004 
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(b) 

(a) 

FIG. 2. Rapidity gaps in deep inelastic scattering: (a) due to soft colors effects; (b) due to 
pomeron exchange. 
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jet 

(b) 

(a) 

FIG. 3. Rapidity gap between jets: (a) due to Pomeron exchange; (b) due to soft color inter- 
actions. 
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FIG. 4. Theoretical fit to the D0 preliminary data Notice that each point shown has a 
different Et bin, with the data being plotted at the mean value of Et for each bin, which is given 
in Table | 
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FIG. 6. Same results of Fig. ^ as function of the adimensional variable xt = Et/^/s- 
comparison we exhibit the recent D0 data [ffl]. 
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FIG. 7. Comparison with the preliminary D0 data |l| of our model prediction for the rapidity 
gap fraction as function of the size of the gap. 
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